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The effect of particle-hole fluctuations for the BCS-BEC crossover is investigated by use of func- 
tional renormalization. We compute the critical temperature for the whole range in the scattering 
length a. On the BCS side for small negative a we recover the Gorkov approximation, while on 
the BEC side of small positive a the particle-hole fluctuations play no important role, and we find 
a system of interacting bosons. In the unitarity limit of infinite scattering length our quantitative 
estimate yields T c /Tp = 0.264. We also investigate the crossover from broad to narrow Feshbach 
resonances - for the later we obtain T c /Tf ~ 0.204 for a -1 = 0. A key ingredient for our treatment 
is the computation of the momentum dependent four-fermion vertex and its bosonization in terms 
of an effective bound-state exchange. 



I. INTRODUCTION 

Ultracold gases of two-component fermions near 
a Feshbach resonance show a smooth crossover be- 
tween Bardeen-Cooper-Schrieffer (BCS) superfluidity 
and Bose-Einstein condensation (BEC) of molecules 
The quantitatively precise understanding of BCS-BEC 
crossover physics is a challenge for theory. Exper- 
imental breakthroughs as the realization of molecule 
condensates and the subsequent crossover to a BCS- 
like state of weakly attractively interacting fermions 
have been achieved (2j. Future experimental precision 
measurements could provide a testing ground for non- 
perturbative methods. An attempt in this direction are 
the recently published measurements of the critical tem- 
perature Q and collective dynamics Q. 

A wide range of qualitative features of the BCS-BEC 
crossover is already well described by extended mean- 
field theories which account for the contribution of both 
fermionic and bosonic degrees of freedom @, Q • In the 
limit of narrow Feshbach resonances mean-field theory 
becomes exact 0, H|. Around this limit perturbative 
methods for small Yukawa couplings 0] can be applied. 
Using e-expansion d or 1/JV -expansion [l(| techniques 
one can go beyond the case of small Yukawa couplings. 

Quantitative understanding of the crossover at and 
near the resonance has been developed through numerical 
calculations u sing various quantum Monte-Carlo (QMC) 
methods [ill Il2l . [T3L H3 |. Computations of the com- 
plete phase diagram have been performed from functional 
field-theoretical techniques, in particular from f-matrix 
approaches [ll|, Dyson-Schwinger equations 0, 2- 
Partice Irreducible methods 1 171 1 , a nd renormalization- 
group flow equations [H, [H, [2(J [2l| . These unified pic- 
tures of the whole phase diagram^], [l(| EH EE ES HH 
20. 21], however, do not yet reach a similar quantitative 
precision as the QMC calculations. 

In this paper we discuss mainly the limit of broad Fesh- 
bach resonances for which all thermodynamic quantities 
can be expressed in terms of two dimensionless parame- 
ters, namely the temperature in units of the Fermi tem- 



perature T /Tp and the concentration c = akp. Here, a 
is the scattering length and the density of atoms is used 
to define kp and Tp via n — k F /(3ir 2 ) and Tp = k 2 F . 
(We use natural units with h = ks = 2M = 1.) In the 
broad resonance regime, macroscopic observables are to 
a large extent independent of the concrete microscopic 
physical realization, a property referred to as universal- 
ity 0, Hot fl9j . This universality includes the unitarity 
regime where the scattering length diverges, a -1 = 
[22| . however it is not restricted to that region. 

For small and negative scattering length c" 1 < 0, |c| <C 
1 (BCS side), the system can be treated with perturba- 
tive methods. However, there is a significant decrease 
in the critical temperature as compared to the original 
BCS result. This was first recognized by Gorkov and 
Melik-Barkhudarov 23]. The reason for this correction 
is a screening effect of particle-hole fluctuations in the 
medium (24|. There has been no systematic analysis of 
this effect in approaches encompassing the full BCS-BEC 
crossover so far. 

In this paper, we present an approach using an exact 
flow equation for the "average action" [25[ (or "flowing 
action"). We include the effect of particle-hole fluctua- 
tions and recover the Gorkov correction on the BCS side. 
We calculate the critical temperature for the second- 
order phase transition between the normal and the su- 
pcrfluid phase throughout the whole crossover. 

We also calculate the critical temperature at the point 
a -1 = for different resonance widths AB. As a function 
of the microscopic Yukawa coupling h\, we find a smooth 
crossover between the exact narrow resonance limit and 
the broad resonance result. The resonance width is con- 
nected to the Yukawa coupling via AB = h\/(87T^MO,b) 
where is the magnetic moment of the bosonic bound 
state and ab is the background scattering length. 

The paper is organized as follows. In Sect. |TT] we 
explain our microscopic model and the flow equation 
method we use to investigate it. Section HTT1 reviews the 
effect of particle-hole fluctuations on a BCS-like super- 
fluid. Section [TVl generalizes this to the whole BCS-BEC 
crossover and explains how we incorporate particle-hole 
fluctuations in our renormalization-group treatment us- 
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ing bosonization. We present our results for the critical 
temperature and the phase diagram in Sees. fVl and [VTl 
Finally, Sect. IVIII describes the crossover from broad to 
narrow resonances and we draw conclusions in Sect. IVIIII 



II. MODEL AND METHOD 

A. Microscopic model 

We investigate the functional integral for the BCS- 
BEC crossover of an ultracold gas of fermionic atoms 
near a Feshbach resonance. We start with a micro- 
scopic action including a two-component Grassmann field 
ip = (ipi, "02) i describing fermions in two hyperfine states. 
Additionally, we introduce a complex scalar field <f> as the 
bosonic degrees of freedom. In different regimes of the 
crossover, it can be seen as a field describing molecules, 
Cooper pairs or simply an auxiliary field. Using the re- 
sulting two-channel model we can describe both narrow 
and broad Feshbach resonances in a unified setting. Ex- 
plicitly, the microscopic action at the ultraviolet scale A 
reads 

ShM = J ' dr J d 3 x{^(d T - A-^ip 

-h A ((fih.ih + h.a)} , (1) 

where we choose nonrelativistic natural units with h = 
ks = = 1, with M the mass of the atoms. The sys- 
tem is assumed to be in thermal equilibrium, which we 
describe using the Matsubara formalism. In addition to 
the position variable x, the fields depend on the imag- 
inary time variable r which parameterizes a torus with 
circumference 1/T. The variable [i is the chemical po- 
tential. The Yukawa coupling h couples the fermionic 
and bosonic fields. It is directly related to the width of 
the Feshbach resonance. The parameter v depends on 
the magnetic field and determines the detuning from the 
Feshbach resonance. Both h and v get renormalized by 
fluctuations, and the microscopic values h\, and v\ have 
to be determined by the properties of two body scattering 
in vacuum. For details, we refer to Qji |26|. 

More formally, the bosonic field <j> appears quadrati- 
cally in the microscopic action in Eq. (TTJ) . The functional 
integral over <fi can be carried out. This shows that our 
model is equivalent to a purely fermionic theory with an 
interaction term 

<Sint = / (- 5 -, H : rWiVJV'lOl) 

x^;(p' 2 )MP2) 5(pi +P2-p' 1 -p' 2 ), (2) 
where p = (j>o,p) and the classical inverse boson propa- 



gator is given by 

P<t>(<l) = *<?o + y + - 2/z . (3) 

On the microscopic level the interaction between the 
fermions is described by the tree level expression 

h 2 

A^.eflF = i — . (4) 

— u) + 2<r + va — 2// 

Here, u> is the real-time frequency of the exchanged boson 
(j). It is connected to the Matsubara frequency qg via 
analytic continuation u — —iqo- Similarly, q = pi+P2 is 
the center of mass momentum of the scattering fermions 
ipi and ip2 with momenta p\ and pi , respectively. 

In this paper we will mainly discuss the limit of broad 
Feshbach resonances, which is realized in current exper- 
iments, e.g. with 6 Li and 40 K. This corresponds to the 
limit h — * oo, for which the microscopic interaction be- 
comes pointlike, with strength —h 2 /v A . For broad Fesh- 
bach resonances, one has a far going universality. Macro- 
scopic quantities are independent of the microscopic de- 
tails and can be expressed in terms of only a few pa- 
rameters. In our case this is the two-body scattering 
length a or, at finite density, the concentration c = a/cp, 
where the Fermi momentum is related to the density by 
kp = (Sn 2 !^) 1 / 3 . At nonzero temperature, an additional 
parameter is given by T/Tp, where Tp is the Fermi tem- 
perature. 

B. Flow equations for the effective action 

The functional renormalization group connects the mi- 
crophysics, as introduced in Sect III Al to macrophysics 
and therefore to observable thermodynamics, by means 
of a non-perturbative flow equation. For this description, 
we start with a functional integral representation of the 
grand canonical partition function 

Z = J Vxe~ sm . (5) 

Here, x collects bosonic and fermionic degrees of freedom. 
In our case it can be written as a vector of the form 

X(x) = U>1 (x), 02 0), $1 (x), 02 0), i>t (x), 02 0)) . (6) 

where we have decomposed the complex bosonic field into 
its two real components, 

4>{x) = -=(4> 1 (x)+i4> 2 (x)). (7) 

We generalize equation (O by introducing a source term 
J J\ and an infrared cutoff term AS/- [%] , so that we ob- 
tain 
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Z k [J] = 



- P w h [j] _ 



= / T>xe 



~ -S[x]-AS k [x]+f Jx 



(8) 



Here, the first equation includes a definition of the in- 
frared regulated functional Wk[J]. The cutoff term 
ASfc[x] is quadratic in the fields. In Fourier space it 
reads 



&s k [x] = \j x f (q)Rk(q)x(q) 



(9) 



where the infrared cutoff function R k (q) has the proper- 
ties 



Rk(q) 
Rk(q) 

Rk(q) 



oo for k — > A , 

e for M^ , 
k 

for — -> . 

\q\ 



(10) 



For k — the functional Zk[J) generates the n-point 
functions, while VFfc[J] is the generating functional for 
the connected correlation functions. The average action 
is defined as a modified Legendre transform of W k [J] , 



coupled differential equations for running couplings. In 
perturbation theory, the solution of Eq. (fT3"|) corre- 
sponds to the computation of infinitely many Feynman 
diagrams. Finding exact solutions to Eq. p3[) for non- 
trivial theories is not possible in practice, but one can 
use truncations in the space of possible functionals to find 
approximate solutions. Such approximations do not have 
to rely on the existence of a small expansion parameter 
such as the interaction strength and they are therefore of 
a non-perturbative nature. For reviews of the functional 
renormalization group method see [13, [H, ■ 

From the full effective action r[x] = Ffe = o[x] one can 
derive all macroscopic properties of the system under 
consideration. For example, the thermodynamic proper- 
ties can be obtained from the grand canonical partition 
function Z or the corresponding grand canonical poten- 
tial <I>g = —T\nZ. It is related to the effective action 
via 



Y[ Xeq ] = $c/T. (15) 
Here \eq is the solution of the field equation 



r[x] 



= o. 



(16) 



-w k [j] 



Jj X ^~AS k [x}, (11) 



with x = (x)J = SWk/SJ. The average action has the 
important property that it interpolates between the mi- 
croscopic action for k — A and the quantum effective 
action, when k = and the cutoff is absent. For k ^ 
the average action Tk[x] is the coarse grained free energy 



S[x] (k 
Tlx] (k 



A) 
0). 



(12) 



The scale dependence of the average action is given by 
the flow equation [25| 



d k T k [ X } = ^STr 



.(2) 



R k d k Rk 



(13) 



Here, the STr operation involves an integration over mo- 
menta as well as a summation over internal indices. We 
employ the second functional derivative of 



(rfM), 



■r fc [x] 



(14) 



Equation (|13p is an exact equation and the starting point 
of our investigations. It is a differential equation for a 
functional which translates to a system of infinitely many 



Since the effective action is the generating functional of 
the one-particle-irreducible (1PI) correlation functions, 
we can also derive dynamical properties from T[x\. For 
example, for a homogeneous field x cq , the full propagator 
G(q) is obtained from 



(T^) ij ( Pu p 2 ) = (G- 1 ) l ,( Pl )6( Pl -p 2 ). 



(17) 



In this work, we solve the flow equation (|13[) approx- 
imately by using a truncation in the space of possible 
functionals T k . More explicitly, our truncation reads 



l/T 



dr I d 3 x \ v> f (a 



+ 4>*(Z <p d T --A^A)0 + U(p, f i) 



(18) 



Here the effective potential U(p,fi) contains no deriva- 
tives and is a function of p = tf>*<f> and p. Besides 
the couplings parameterizing U (see below) our trunca- 
tion contains three further /c-dependent ( "running" ) cou- 
plings A<p, Z^ and h. The truncation in Eq. (fT8)) can 
be motivated by a systematic derivative expansion and 
analysis of symmetry constraints (Ward identities), see 
[Hii, [H. The truncation in Eq. does not yet 
incorporate the effects of particle-hole fluctuations and 
we will come back to this issue in Sect. IIV1 Eq. 



In terms of renormalized fields 
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renormalized couplings Z$ — Z^/A^, h — hj ' \J A^ and 
effective potential U(p,p) — U(p 7 p), Eq. (JTBJ) reads 

r fe [ X ] = J ' drj d 3 x^(d T -A-p)i> 

+ ^{Z^d T -\^)<P + U{p^) 

- M^Vi^+^i*)}- (19) 

For the effective potential, we use an expansion around 
the fc-dependent location of the minimum po(k) and the 
fc-independent value of the chemical potential po that 
corresponds to the physical particle number density n. 
We determine pa(k) and p by the requirements 

(d p U)(p (k),p ) = for all k 
-(d fl U){p ,p )=n atfc = 0. (20) 



reads in terms of the bare fields 

AS k = J |^t(p) ( sign (^ _ A1 )fc 2 - {f - p)) 
x 9 (k 2 - \f - p\) ifr(p) 
+4>*(p)A^ (k 2 ~f/2) 9 (k 2 -f/2) ftp)}. 

(22) 

For the fermions it regularizes fluctuations around the 
Fermi surface, while for the bosons fluctuations with 
small momenta are suppressed. The choice of AS k in 
Eq. (|22|) is an optimized choice in the sense of [1^, [3l[ . 

The flow equations without the effect of the particle- 
hole exchange have been discussed extensively elsewhere 
(l9l I20I ]. Even though we use here a different cutoff 
and a different treatment of the density, we will not give 
them explicitly and rather concentrate on the effects of 
particle-hole exchange diagrams. 



More explicitly, we employ a truncation for U (p, p) of the 
form 

U(p,p) = U(p ,po) -n k (p- p ) 

+ (m 2 + a(p - p ))(p - p ) 

+ \X(P- P0 ) 2 - (21) 

In the symmetric or normal gas phase, we have po = 0, 
while in the regime with spontaneous symmetry break- 
ing, we have m 2 = 0. The atom density n = —dU/dp 
corresponds to n k in the limit k — > 0. 

In total, we have the running couplings m 2 (fc), A(fc), 
a(k), nk, Zrf,{k) and h(k). (In the phase with sponta- 
neous symmetry breaking m 2 is replaced by po-) In ad- 
dition, we need the anomalous dimension r\ = — kdk^A^. 
We project the flow of the average action T k on the 
flow of these couplings by taking appropriate (functional) 
derivatives on both sides of Eq. (fT5|) . We thereby obtain 
a set of coupled nonlinear differential equations which 
can be solved numerically. 

At the microscopic scale k = A the initial values of 
our couplings are determined from Eq. (|12p . This gives 
m 2 (A) = v A - 2p , po = 0, A (A) = 0, Z^A) = 1, 
h(A) = h\, a(A) = —2 and ija = 3TT 2 p n 8(p n ). The initial 
values va and h\ can be connected to the two particle 
scattering in vacuum close to a Feshbach resonance. For 
this purpose one follows the flow of m 2 (k) and h(k) in 
vacuum, i.e. po = T = n = and extracts the renor- 
malized parameters m 2 — m 2 (k = 0), h = h(k = 0). 
The scattering length a obeys a = —h 2 /(8Trm 2 ) and the 
renormalized Yukawa coupling h determines the width of 
the resonance as discussed in Sect. IVIIl Broad Feshbach 
resonances with large h become independent of h. 

The infrared cutoff we use is purely space-like and 



III. PARTICLE-HOLE FLUCTUATIONS 

The BCS theory of superfluidity in a Fermi gas of 
atoms is valid for a small attractive interaction between 
the fermions [32l |. In a renormalization group setting, 
the features of BCS theory can be described in a purely 
fcrmionic language. The only scale dependent object is 
the fermion interaction vertex X^. The flow depends on 
the temperature and the chemical potential. For positive 
chemical potential (p > 0) and small temperatures T, 
the appearance of pairing is indicated by the divergence 
of 

In general, the interaction vertex is momentum depen- 
dent and represented by a term 

I\, = / \^{p' ll pi,p' 2 ,p 2 ) 

J Pl,P2,Pi,P'2 

xrM)MpiW 2 (P2)MP2) (23) 

in the effective action. In a homogeneous situation, mo- 
mentum conservation restricts the expression in Eq. (|23[) 
to three independent momenta, A^, ~ 5(p' 1 +p 2 ~Pi —P2)- 
The flow of Aip has two contributions which are depicted 
in Fig. [TJ The first diagram describes particle-particle 
fluctuations. For p > its effect increases as the temper- 
ature T is lowered. For small temperatures T < T c .bcs 
the logarithmic divergence leads to the appearance of 
pairing, as A^, — > 00. 

In the purely fermionic formulation the flow equation 
for A,/, has the general form [28|, [33| 

d k \% = A^XK . ( 24 ) 

with a, /9, 7 denoting momentum as well as spin labels. 
A numerical solution of this equation is rather involved 
due to the rich momentum structure. The case of the 
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FIG. 1: Running of the momentum dependent vertex Ay,. 
Here dk indicates derivatives with respect to the cutoff terms 
in the propagators and does not act on the vertices in the 
depicted diagrams. We will refer to the first loop as the 
particle-particle loop (pp-loop) and to the second one as the 
particle-hole loop (ph-loop). 



attractive Hubbard model in two dimensions, which is 
close to our problem, has recently been discussed in [34j . 
The BCS approach concentrates on the pointlike cou- 
pling, evaluated by setting all momenta to zero. For 
k — > 0, fio — > 0, T — > and n — > this coupling is 
related the scattering length, a = -^\^{pi = 0). In the 
BCS approximation only the first diagram in Fig. [T] is 
kept, and the momentum dependence of the couplings 
on the right-hand side of Eq. (|2"4")l is neglected, by re- 
placing by the pointlike coupling evaluated at zero 
momentum. In terms of the scattering length a, Fermi 
momentum kp and Fermi temperature Tp, the critical 
temperature is found to be 



To 



a61e V(2afc F ) 



(25) 



This is the result of the original BCS theory. However, 
it is obtained by entirely neglecting the second loop in 
Fig. [TJ which describes particle-hole fluctuations. At 
zero temperature the expression for this second diagram 
vanishes if it is evaluated for vanishing external momenta. 
Indeed, the two poles of the frequency integration are 
always either in the upper or lower half of the complex 
plane and the contour of the frequency integration can 
be closed in the half plane without poles. 

The dominant part of the scattering in a fermion gas 
occurs, however, for momenta on the Fermi surface rather 
than for zero momentum. For non-zero momenta of the 
" external particles" the second diagram in Fig. [1] - the 
particle-hole channel - makes an important contribution. 

Setting the external frequencies to zero, we find that 
the inverse propagators in the particle-hole loop arc 



P^(q) = m + {q-pif -m, 



(26) 



planes so that there is a nonzero contribution even for 
T = 0. 

To include the effect of particle-hole fluctuations one 
could try to take the full momentum dependence of the 
vertex A^, into account. However, this leads to compli- 
cated expressions which are hard to solve even numeri- 
cally. One therefore often restricts the flow to the running 
of a single coupling A,/, by choosing an appropriate projec- 
tion prescription to determine the flow equation. In the 
purely fermionic description with a single running cou- 
pling \tp, this flow equation has a simple structure. The 
solution for A^ 1 can be written as a contribution from 
the particle-particle (first diagram in Fig. [TJ pp-loop) 
and the particle-hole (second diagram, ph-loop) channels 



1 



1 



A^(fc = 0) A^(fe = A) 



pp-loop + ph-loop . (28) 



Since the ph-loop depends only weakly on the tempera- 
ture, one can evaluate it at T — and add it to the initial 
value \^{k — A) -1 . Since T c depends exponentially on 
the "effective microscopic coupling" 



(a: 



eff \ 

0,aJ 



= X,^(k = A) 1 + ph-loop , 



(29) 



any shift in ( A J results in a multiplicative factor 

for T c . The numerical value of the ph-loop and therefore 
of the correction factor for T c /Tp depends on the precise 
projection description. 

Let us now choose the appropriate momentum config- 
uration. For the formation of Cooper pairs, the relevant 
momenta lie on the Fermi surface, 



-Q -9 -»/? -*F2 
Pl=P 2 = Pi = P2 = A* , 



and point in opposite directions 



Pi 



"P2, Pi = ~P 2 ■ 



(30) 



(31) 



This still leaves the angle between p\ and p{ unspeci- 
fied. Gorkov's approximation uses Eqs. ([30)) and (|3Tj) 
and projects on the s-wave by averaging over the an- 
gle between p\ and p[. One can shift the loop momen- 
tum such that the internal propagators depend on q 2 and 
(q + Pi — Pi) 2 - In terms of spherical coordinates the first 
propagator depends only on the magnitude of the loop 
momentum q 2 = (f 1 , while the second depends addition- 



ally on the transfer momentum p 2 
angle a between q and (pi — p[), 



3 (pl - Pi) 2 and the 



and 



pAq) = i( lo + (3-P2) - P 



(27) 



Depending on the value of the momenta p\ and p^ there 
are now values of the loop momentum q for which the 
poles of the frequency integration are in different half 



(q+ Pi - P[f = I 2 + 4p 2 + 4gpcos(a) . (32) 

Performing the loop integration involves the integration 
over q 2 and the angle a. The averaging over the angle be- 
tween pi and p[ translates to an averaging over p 2 . Both 
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can be done analytically [2J| for the fermionic particle- 
hole diagram and the result gives the well-known Gorkov 
correction to BCS theory, resulting in 



1 



( 4e )i/3 



T, 



c,BCS 



1 

%2 



c.BCS ■ 



(33) 



In this paper we will use a numerically simpler pro- 
jection by choosing p[ = pi, and pi — p^ without an 
averaging over the angle between p[ and p\. The size 
of p 2 = p\ is chosen such that the one-loop result repro- 
duces exactly the result of the Gorkov correction, namely 
p = 0.732Q^fjl. Choosing different values of p demon- 
strates the dependence of T c on the projection procedure 
and may be taken as an estimate for the error that arises 
from the limitation to one single coupling A^ instead of 
a momentum dependent function. 



IV. BOSONIZATION 

In Sec. |TT] we describe an effective four-fermion in- 
teraction by the exchange of a boson. In this picture 
the phase transition to the superfluid phase is indicated 
by the vanishing of the bosonic "mass term" m 2 = 0. 
Negative m 2 leads to the spontaneous breaking of U(l)- 
symmetry, since the minimum of the effective potential 
occurs for a nonvanishing superfluid density po > 0. 



Pi 



P2 





t 



Pl P2 

Pi P2 

FIG. 2: Flow of the boson propagator. 

For m 2 > we can solve the field equation for the 
boson <j) as a functional of ip an d insert the solution 
into the effective action. This leads to an effective four- 
fermion vertex describing the scattering ipi(pi)ip2{p2) —> 

MpDMp?) 



A 



V>,eff : 



i(pi +P2)o + h(Pi + P2) 2 + 



(34) 



To investigate the breaking of U(l) symmetry and the 
onset of superfluidity, we first consider the flow of 
the bosonic propagator, which is mainly driven by the 
fermionic loop diagram. For the effective four-fermion in- 
teraction this accounts for the particle-particle loop (see 



r.h.s. of Fig. [2]). In the BCS limit of a large microscopic 
m\ the running of m 2 for k — > reproduces the BCS 
result HI]. 

The particle-hole fluctuations are not accounted for by 
the renormalization of the boson propagator. Indeed, we 
have neglected so far that a term 



(35) 



in the effective action is generated by the flow. This holds 
even if the microscopic pointlike interaction is absorbed 
by a Hubbard-Stratonovich transformation into an effec- 
tive boson exchange such that A,/, (A) = 0. The strength 
of the total interaction between fermions 



A,/,.eff = 



+ A 



(36) 



i{Pi +P2)q + 5 (Pi +P2) Z + m l 
adds A,/, to the piece generated by boson exchange. In 

P2 P2 



fie 



Pi 



Pi 



FIG. 3: Box diagram for the flow of the four-fermion interac- 
tion. 

the partially bosonized formulation, the flow of A^ is gen- 
erated by the box-diagrams depicted in Fig. [31 We may 
interpret these diagrams and establish a direct connec- 
tion to the particle-hole diagrams depicted in Fig. [1] 
on the BCS side of the crossover and in the microscopic 
regime. There the boson gap m 2 is large. In this case, 
the effective fermion interaction in Eq. (|3G[) becomes 
momentum independent. Diagrammatically, this is rep- 
resented by contracting the bosonic propagator. One can 
see, that the box-diagram in Fig. [3] is then equivalent to 
the particle-hole loop investigated in Sec. IIIII with the 
pointlike approximation A^ je ff — ► — ^ for the fermion 
interaction vertex. As mentioned above, these contribu- 
tions vanish for T = 0, /1 < for arbitrary pi. Indeed, 
at zero temperature, the summation over the Matsubara 
frequencies becomes an integral. All the poles of this in- 
tegration are in the upper half of the complex plane and 
the integration contour can be closed in the lower half 
plane. We will evaluate fifcA^, for pi = p{ = —P2 — —P2, 
\pi\ = p = 0.7326 v / ^(, as discussed in the Sec. IIIII For 
/i > this yields a nonvanishing flow even for T = 0. 

Another simplification concerns the temperature de- 
pendence. While the contribution of particle-particle di- 
agrams becomes very large for small temperatures, this 
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is not the case for particle- hole diagrams. For nonvan- 
ishing density and small temperatures, the large effect 
of particle-particle fluctuations leads to the spontaneous 
breaking of the U(l) symmetry and the associated super- 
fluidity. In contrast, the particle- hole fluctuations lead 
only to quantitative corrections and depend only weakly 
on temperature. This can be checked explicitly in the 
pointlike approximation, and holds not only in the BCS 
regime where Tj\i <C 1, but also for moderate T//z as re- 
alized at the critical temperature in the unitary regime. 
We can therefore evaluate the box-diagrams in Fig. [T]for 
zero temperature. We note that an implicit temperature 
dependence, resulting from the couplings parameterizing 
the boson propagator, is taken into account. 

After these preliminaries, we can now incorporate the 
effect of particle-hole fluctuations in the renormalization 
group flow. A first idea might be to include the additional 
term (|35p in the truncation and to study the effects of 
X^p on the remaining flow equations. On the initial or 
microscopic scale one would have A,/, = 0, but it would 
then be generated by the flow. This procedure, however, 
has several shortcomings. First, the appearance of a local 
condensate would now be indicated by the divergence of 
the effective four-fermion interaction 



In consequence, the flow equations in the symmetric 
regime get modified 



d k h 



dkh\^ - P r p(q)d k v. 
<9fcA</,| t - 2hd k v. 



(40) 
(41) 



Here q is the center of mass momentum of the scattering 
fermions. In the notation of Eq. (|34[) we have q = p\ +P2 
and we will take q = 0, and go — 0- The first term on 
the right hand side in Eq. (|4"Tj) gives the contribution of 
the flow equation which is valid for fixed field <p k . The 
second term comes from the explicit scale dependence of 
cf>k- The inverse propagator of the complex boson field (f> 
is denoted by P<p(q) = A^P^q) = A<t,(m 2 +iZ<f,qo+(f /2), 
cf. Eq. ©. 

We can choose d k v such that the flow of the coupling 
A^, vanishes, i.e. that we have A^ = on all scales. This 
modifies the flow equation for the renormalized Yukawa 
coupling according to 



d k h = d k h\- 4>k -—d k \ ij \^ 



(42) 



A 



V>,eff 



+ A 



(37) 



This might lead to numerical instabilities for large or di- 
verging A^. The simple picture that the divergence of 
A^, c ff is connected to the onset of a nonvanishing expec- 
tation value for the bosonic field (f>o , at least on interme- 
diate scales, would not hold anymore. Furthermore, the 
dependence of the box-diagrams on the center of mass 
momentum would be neglected completely by this proce- 
dure. Close to the resonance the momentum dependence 
of the effective four-fermion interaction in the bosonized 
language as in Eq. (|36|) is crucial, and this might also be 
the case for the particle-hole contribution. 

Another, much more elegant way to incorporate the 
effect of particle-hole fluctuations is provided by the 
method of bosonization [5^, [35|. For this purpose, we use 
scale dependent fields in the average action. The scale 
dependence of Tfc[xfe] is modified b y a term reflecting the 
fc-dependence of the argument \k [35| 



d k T k [x k ] = I ^f^dkXk + ^STr 
0Xk 2 



d k R k 



( 38 ) 

For our purpose it is sufficient to work with scale de- 
pendent bosonic fields <fi and keep the fermionic field ip 
scale independent. In practice, we employ bosonic fields 
and 4> k with an explicit scale dependence which reads 
in momentum space 



d k (t) k (q) 
dk^Uq) 



(39) 



with di h 7 the contribution without bosonization and 

® k ^\d> k §i ven by the box diagram in Fig. [3l Since X^ 
remains zero during the flow, the effective four-fermion 
interaction A^, j0 ff is now purely given by the boson ex- 
change. However, the contribution of the particle-hole 
exchange diagrams is incorporated via the second term 
in Eq. glD. 

In the regime with spontaneously broken symmetry we 
use a real basis for the bosonic field 



1 

V2 V 



(43) 



where the expectation value </>o is chosen to be real with- 
out loss of generality. The real fields </>i and </>2 then de- 
scribe the radial and the Goldstone mode, respectively. 
To determine the flow equation of h, we use the projec- 
tion description 



d k h = i\/2Vt 



S 



Sfo(o) (SVi(o) smo) 



d k T k 



(44) 



with the four volume f2 = A L. Since the Goldstone 
mode has vanishing "mass" , the flow of the Yukawa cou- 
pling is not modified by the box diagram (Fig. [3]) in the 
regime with spontaneous symmetry breaking. 

We emphasize that the non-perturbative nature of the 
flow equations for the various couplings provides for a re- 
summation similar to the one in Eq. (|2"5)) , and thus goes 
beyond the treatment by Gorkov and Melik-Barkhudarov 
[23j which includes the particle-hole diagrams only in a 
perturbative way. Furthermore, the inner bosonic lines 
h 2 / P^q) in the box-diagrams represent the center of 
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mass momentum dependence of the four-fermion vertex. 
This center of mass momentum dependence is neglected 
in Gorkov's pointlike treatment, and thus represents a 
further improvement of the classic calculation. Actually, 
this momentum dependence becomes substantial - and 
should not be neglected in a consistent treatment - away 
from the BCS regime where the physics of the bosonic 
bound state sets in. Finally, we note that the trunca- 
tion (fT9]l supplemented with ([35)1 closes the truncation 
to fourth order in the fields except for a fermion-boson 
vertex X^^i/j' t/j(j)* <p which plays a role for the scattering 
physics deep in the BEC regime [36] but is not expected 
to have an important impact on the critical temperature 
in the unitarity and BCS regimes. 

V. CRITICAL TEMPERATURE 

To obtain the flow equations for the running couplings 
of our truncation Eq. (|19[) we use projection prescriptions 
similar to Eq. (|44[) . The resulting system of ordinary cou- 
pled differential equations is then solved numerically for 
different chemical potentials p and temperatures T. For 
temperatures sufficiently small compared to the Fermi 
temperature T F = (37r 2 n) 2 / 3 , T/T F « 1 we find that 
the effective potential U at the macroscopic scale k = 
develops a minimum at a nonzero field value po > 0, 
dpU(po) = 0. The system is then in the superfluid phase. 
For larger temperatures we find that the minimum is at 
po = and that the "mass parameter" m 2 is positive, 
m 2 = d p U(0) > 0. The critical temperature T c of this 
phase transition between the superfluid and the normal 
phase is then defined as the temperature where one has 

po = 0, d p U(0) = at k = 0. (45) 

Throughout the whole crossover the transition po — * 
is continuous as a function of T demonstrating that the 
phase transition is of second order. 

In Fig. Q] we plot our result obtained for the critical 
temperature T c and the Fermi temperature Tp as a func- 
tion of the chemical potential p at the unitarity point 
with a -1 = 0. From dimensional analysis it is clear that 
both dependencies are linear, T c , Tp ~ p, provided that 
non-universal effects involving the ultraviolet cutoff scale 
A can be neglected. That this is indeed found numeri- 
cally can be seen as a nontrivial test of our approximation 
scheme and the numerical procedures as well as the uni- 
versality of the system. Dividing the slope of both lines 
gives T c /Tp = 0.264, a result that will be discussed in 
more detail below. We emphasize that part of the poten- 
tial error in this estimates is due to uncertainties in the 
precise quantitative determination of the density or Tp. 

VI. PHASE DIAGRAM 

The effect of the particle-hole fluctuations shows most 
prominently in the result for the critical temperature. 



T 




0.0 0.2 0.4 0.6 0.8 



FIG. 4: (Color online) Critical temperature T c (boxes) and 
Fermi temperature Tp = (37r 2 n) 2/ ' 3 (triangles) as a function 
of the chemical potential p. For convenience the Fermi tem- 
perature is scaled by a factor 1/5. We also plot the linear fits 
T c — 0.39/i and Tf = lA8p. The units are arbitrary and we 
use A = e 7 . 

With our approach we can compute the critical tempera- 
ture for the phase transition to superfluidity throughout 
the crossover. The results are shown in Fig. [5J We 
plot the critical temperature in units of the Fermi tem- 
perature T c /Tp as a function of the scattering length 
measured in units of the inverse Fermi momentum, i. e. 
the concentration c = akp. We can roughly distinguish 

T 

0.35^ , — r 7 — I ■ — ■ 1 — ^ ■ — n 

0.30 : : 




FIG. 5: (Color online) Dimensionless critical temperature 
T c /Tf as a function of the inverse concentration c _1 = 
(o/cf) -1 . The black solid line includes the effect of particle- 
hole fluctuations. We also show the result obtained when 
particle- hole fluctuations are neglected (dot-dashed line). For 
comparison, we plot the BCS result without (left dotted line) 
and with Gorkov's correction (left dashed). On the BEC side 
with c _1 > 1 we show the critical temperature for a gas of 
free bosonic molecules (horizontal dashed line) and a fit to 
the shift in T c for interacting bosons, AT C ~ c (dotted line on 
the right). The black solid dot gives the QMC results [infill]. 

three different regimes. On the left side, where cT 1 < —1, 
the interaction is weakly attractive. Mean field or BCS 
theory is qualitatively valid here. In Fig. [5] we denote 
the BCS result by the dotted line on the left (c" 1 < 0). 
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However, the BCS approximation has to be corrected by 
the effect of particle-hole fluctuations, which lower the 
value for the critical temperature by a factor of 2.2. This 
is the Gorkov correction (dashed line on the left side in 
Fig. [5]). The second regime is found on the far right 
side, where the interaction again is weak, but now we 
find a bound state of two atoms. In this regime the sys- 
tem exhibits Bose-Einstein condensation of molecules as 
the temperature is decreased. The dashed horizontal line 
on the right side shows the critical temperature of a free 
Bose-Einstein condensate of molecules. In-between there 
is the unitarity regime, where the two-atom scattering 
length diverges (c _1 — > 0) and we deal with a system of 
strongly interacting fcrmions. 

Our result including the particle-hole fluctuations is 
given by the solid line. This may be compared with a 
functional renormalization flow investigation without in- 
cluding particle- hole fluctuations (dot-dashed line) [l9j |. 
For c — > 0_ the solid line of our result matches the BCS 
theory including the correction by Gorkov and Melik- 
Barkhudarov 12311 . 



To 
T F 



7/3 



2 tt/(2c) _ o.28e 7r/(2c) . 



(46) 



In the regime c _1 > —2 we see that the non-perturbative 
result given by our RG analysis deviates from Gorkov's 
result, which is derived in a perturbative setting. 

On the BEC-side for very large and positive c _1 our 
result approaches the critical temperature of a free Bose 
gas where the bosons have twice the mass of the fermions 
Mb = 2M. In our units the critical temperature is then 



T, 



e,BEC 



2tt 



0.218. 



Tf [6tt 2 C(3/2)] 2/3 
For c — > 0+ this value is approached in the form 



T, 



Cc.BEC 



T 



c.BEC 



1/3 

Ka M n M 



a (67T 2 ) 1 / 3 ' 



(47) 



(48) 



Here, tim = n/2 is the density of molecules and clm is 
the scattering length between them. For the ratio ajf/a 
we use our result om/o = 0.718 obtained from solving 
the flow equations in vacuum, i. e. at T = n = 0, see 
also This has to be compared to the result ob- 

tained from solving the corresponding Schrodinger equa- 
tion which gives clm /a = 0.6 [13] • For the coefficients 
determining the shift in T c compared to the free Bose 
gas we find k — 1.55. 

For c _1 > 0.5 the effect of the particle- hole fluctuations 
vanishes. This is expected since the chemical potential 
is now negative fi < and there is no Fermi surface any 
more. Because of that there is no difference between the 
new curve with particle-hole fluctuations (solid in Fig. 
[5]) and the one obtained when particle-hole contributions 
are neglected (dot-dashed in Fig. [5]). Due to the use of 
an optimized cutoff scheme and a different computation 
of the density our results differ slightly from the ones 
obtained in [lj|. 



In the unitary regime (c _1 w 0) the particle-hole fluc- 
tuations still have a quantitative effect. We can give an 
improved estimate for the critical temperature at the res- 
onance (c" 1 = 0) where we find T c /Tp — 0.264. Results 
from quantum Monte Carlo simulations are T c /Tp = 0.15 
[H El and T C /T F = 0.245 [Q. The measurement 
by Luo et al. 3] in an optical trap gives T c /Tp = 
0.29(+0.03/ — 0.02), which is a result based on the study 
of the specific heat of the system. 



VII. CROSSOVER TO NARROW RESONANCES 

Since we use a two channel model (Eq. |T])) we can 
not only describe broad resonances with h\ — > 00 but 
also narrow ones with h\ —* 0. This corresponds to a 
nontrivial limit of the theory which can be treated exactly 
0, H[ • In the limit h\ — > the microscopic action Eq. ((J) 
describes free fermions and bosons. The essential feature 
is, that they are in thermodynamic equilibrium so that 
they have equal chemical potential. (There is a factor 2 
for the bosons since they consist of two fermions.) For 
vanishing Yukawa coupling the theory is Gaussian 
and the macroscopic propagator equals the microscopic 
propagator. There is no normalization of the "mass"- 
tcrm to 2 so that the detuning parameter in Eq. fT} is 
v = ^m(B - Bq) and 



= fM M (B-B )-2fM. 



(49) 



To determine the critical density for fixed temperature, 
we have to adjust the chemical potential \i such that the 
bosons are just at the border to the superfluid phase. For 
free bosons this implies m 2 = and thus 

H= l pH[ (B-B ) = -^-h 2 a-\ (50) 

2 16-7T 

In the last equation we use the relation between the de- 
tuning and the scattering length 



(51) 



8TTfJ, M (B - Bq 



The critical temperature T c is now determined from the 
implicit equation 



d 3 p 



-J— ff 2 -, 



n. (52) 



While the BCS-BEC crossover can be studied as a 
function of B — Bq or /1, Eq. (|5l"|) implies that for h\ — > 
a finite scattering length a requires B — > Bo- For all 
c ^ the narrow resonance limit implies for the phase 

transition B = Bq and therefore fi = 0. (A different 

concentration variable c mec j was used in 0, ITa | , such that 
the crossover could be studied as a function of c mc( j in 
the narrow resonance limit, see the discussion at the end 
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of this section.) For /i = Eq. (|52j) can be solved ana- 
lytically and gives 



To 
T F 



Ay/2 



2/3 



3(3 + \/2)7ri/2 C (3/2) 



0.204. 



(53) 



This result is confirmed numerically by solving the flow 
equations for different microscopic Yukawa couplings h\ 
and taking the limit h\ — > 0. In Fig. [6l we show the crit- 
ical temperature T c /Tf as a function of the dimension- 
less Yukawa coupling h^/y/kp in the "unitarity limit" 
c _1 = (solid line). For small values of the Yukawa cou- 
pling, ht\j \fkp < 2 we enter the regime of the narrow 
resonance limit and the critical temperature is indepen- 
dent of the precise value of h\. The numerical value 
matches the analytical result T c /Tp w 0.204 (dotted line 
in Fig. [6]). For large Yukawa couplings, h\/y/kp > 40, 
we recover the result of the broad resonance limit as ex- 
pected. In between there is a smooth crossover of the 
critical temperature from narrow to broad resonances. 

T c /Tf 



0.26 



0.24 



0.22- 



0.20- 



0.1 



10 



100 



1000 



Ha/ Vkp 



FIG. 6: (Color online) The critical temperature divided by 
the Fermi temperature T c /Tf as a function of the dimension- 
less Yukawa coupling h,A/Vk~F for c _1 = (solid line). One 
can clearly see the plateaus in the narrow resonance limit 
(Tc/Tf ~ 0.204, dotted line) and in the broad resonance limit 
(T c /Tf « 0.264, dashed line). 

We use here a definition of the concentration c = ak F 
in terms of the vacuum scattering length a. This has the 
advantage of a straightforward comparison with experi- 
ment since a -1 is directly related to the detuning of the 
magnetic field B — B , and the "unitarity limit" cr 1 = 
precisely corresponds to the peak of the resonance B — 
Bq. However, for a nonvanishing density other definitions 
of the concentration parameter are possible, since the ef- 
fective fermion interaction \-^, e g depends on the density. 
For example, one could define for n / a "in medium 
scattering length" a = A^ c ff/(87r), with \ip. e $ = —h 2 /m 2 
evaluated for T — but n ^ @. The corresponding 
"in medium concentration" c mcc \ = ak F would differ from 
our definition by a term involving the chemical potential, 
resulting in a shift of the location of the unitarity limit 
if the latter is defined as c m J d = 0. While for broad res- 
onances both definitions effectively coincide, for narrow 
resonances a precise statement how the concentration is 



defined is mandatory when aiming for a precision com- 
parison with experiment and numerical simulations for 
quantities as T c /Tp at the unitarity limit. For exam- 
ple, defining the unitarity limit by c~ d = would shift 
the critical temperature in the narrow resonance limit to 
T C /T F = 0.185 



VIII. CONCLUSION 

This paper aims at a quantitatively reliable computa- 
tion of the phase diagram for the transition to superfluid- 
ity for non-relativistic fermions. We obtain an estimate of 
the critical temperature T c for the whole range of scatter- 
ing length in the BCS-BEC crossover. A quantitatively 
precise determination of T c for the BCS-BEC crossover is 
a challenge both for theory and experiment. A successful 
match would extend the success of the "universal critical 
physics" to quantities that are "nonuniversal" in the lan- 
guage of critical phenomena as the critical temperature 
itself. On the theoretical side this requires a full control 
of the mapping from microphysics to macrophysics. It 
constitutes a decisive test for non-perturbative methods 
in complex many-body theories. New extended notions of 
universality become visible, which also extend to physics 
away from the critical temperature, as the low tempera- 
ture limit T — > 0. 

For non-relativistic fermions the task of a precise com- 
putation of T c is difficult even for small scattering length. 
The basic reason is that the phase transition itself is 
a non-perturbative phenomenon, linked to an effective 
fermion interaction A,/, )C ff growing to large values. In or- 
der to locate T c precisely, one has to follow the scale- 
dependence of A^eff with sufficient precision. The prob- 
lem arises from the substantial momentum dependence of 
the fluctuation contributions to the effective four-fermion 
vertex. In principle, a consistent one loop approximation 
to the flow should follow the full momentum dependence 
of the vertex according to Eq. (|24|) . With the general 
form T C /T F = Fe 7r ^ 2kFa ^> for ak F -> 0_ this will influ- 
ence the value of the prefactor F. Different projections 
of the flow in Eq. ([24| onto a single coupling A^, , as the 
BCS approximation or the Gorkov projection prescrip- 
tion, result in different approximations for F. 

In this paper we use a projection that coincides with 
the Gorkov approximation in the BCS-limit ak F — > 0_. 
Our method of non-perturbative flow equations allows 
us to extend Gorkov's result to the whole range of a, 
including the "unitarity limit" a — > oo. Already at the 
present stage of truncation our result T c /Tf = 0.264 for 
a — ► oo agrees well with experiment and competes with 
quantum Monte-Carlo simulations. Further extensions of 
the truncation should refine the estimate of T c and also 
permit an indication of the theoretical uncertainty of this 
estimate. 
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